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Abstract 

We present a topological method for the efBcient computer assisted 
verification of the existence of the homoclinic tangency which unfolds 
generically in a one-parameter family of planar maps. The method has 
been applied to the Henon map and the forced damped pendulum ODE. 
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1 Introduction 

The goal of our paper is to describe a method for the verification of the ex- 
istence of the quadratic homoclinic tangency which unfolds generically in a 
one-parameter family of planar maps. This is an important problem problem in 
dynamics, because establishing the existence of the generic homoclinic tangency 
has very deep dynamical consequences, see for example |PT1 IWYll IWY21 IRCj 
and references given there. 

Our paper was inspired by the works of Aral and Mischaikow [XI lAMj , who 
combined some tools from the Conley index theory, the tools of computational 
homology from the CHomP project [CHOMP] and the set oriented numerical 
methods from the GAIO project jGAIOj . into a method for computer assisted 
proof of the existence of the generic homoclinic tangency. Using it they proved in 
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|AM| the existence of the generic homochnic tangency for the dissipative Henon 
map Ha.b{x, y) = {a — x'^ + by, x) for parameter values close to a = 1.4, b = 0.3 
and a — 1.3,6 = —0.3. Their method contains essentially two separate parts: 
first using the Conley index approach they prove the existence of the homoclinic 
tangency for some parameter value (steps 1 to 5 in the terminology used in [AMj ) 
and in the second part (step 6 in |AMj ) they verify some transversality-type 
condition, which implies the genericity of the homoclinic tangency established 
in the first part. The computation times reported in jAMj are around 260 
minutes and 100 minutes on a PowerMac G5 (2GHz) for b = 0.3 and b — —0.3, 
respectively. In these computations the second part took 61 minutes and 24 
minutes, respectively. From these computation times it is quite clear that there 
is a little hope to apply successfully this method to ODEs. 

Our method, also topological and geometric in sprit, is based on the observa- 
tion that the computations done by Arai and Mischaikow in |AMj in the second 
part of their approach should be in principle sufficient to obtain the whole result, 
both the existence of the homoclinic tangency and its genericity. Obviously for 
this end, one should use different tools, ours are essentially those of differential 
topology, but developed earlier in the context of topological dynamics, like the 
covering relations combined with the cone conditions, see [KWZ[ As a test 
case we give a computer assisted proof of the existence of the generic homoclinic 
tangency for the dissipative Henon map with 6 = —0.3 and a 1.3145 in the 
computation time 0.2 sec on the Intel Xeon 5160, 3GIIz processor. This should 
be contrasted with the fact that in [AMj this case took around 100 minutes. 

An application of our method to an ODE requires an efficient rigorous C^- 
solver for ODEs. By this we mean an algorithm for rigorous integration of 
ODEs together with their variational equations up to the second order. Such 
an algorithm has been recently developed by the authors in [WZ2| and is now 
a part of the CAPD library [CAPDj . Using this algorithm we were able to 
prove the existence of the generic homoclinic tangency for the 27r-shift along 
the trajectory of the periodically forced pendulum equation 

X + (3x + sin{x) — cos{t) . (1) 

for f3 w 0.2471 (see Theorem [T51 in Section [7]). The computation time for this 
proof is 30 sec on the Intel Xeon 5160, 3GHz processor. 

It should be mentioned that similar results for the Henon map have been 
obtained by a complex analytic method of Fornaess and Gavosto [FG1|, IFG2| . 
Compared to their method, which depends on the analyticity of maps, our 
method and that of Arai and Mischaikow are rather geometric and topological, 
and are designed so that they can be applied to a wider class of maps. Essen- 
tially, we require a continuous family of diffeomorphisms for which we can 
compute the image of the maps using interval arithmetic. 

The content of the paper may be briefly described as follows: in Sections [2] 
and [3] following mainly [AM] we give basic definitions and restate the problem of 
the existence of the quadratic homoclinic tangency unfolding generically as the 
transversality question for a dynamical system induced by the given one on the 
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projective bundle. In Section S] we discuss how this trans versahty problem can 
be solved using the covering relations linked with the cone conditions. Moreover, 
we derive computable estimates for the dependence on a parameter of (un)stable 
manifolds of the hyperbolic fixed point, which will be later used in the computer 
assisted proofs. In Section Owe illustrate our approach on a toy example, with 
the intention that the reader may see and appreciate some details of the method, 
which are later hardly visible when we report on the computer assisted proof 
for the Henon map and for the forced damped pendulum in Sections [H] and [71 
respectively. 

2 Homoclinic tangency and the projectivization 

2.1 Invariant manifold - basic notation 

Definition 1 Consider the map f : X Z) dom (/) X . 

Let X e X. Any sequence \xk\k^i, where I d 1^ is a set containing and 
for any h < I2 < h in ^ if h,h G I , then I2 G /, such that 

Xo = x, f{xi) = Xi+i, fori,i + leI 

will be called an orbit through x. If I = then we will say that is a 

full backward orbit through x. 

Definition 2 Let X be a topological space and let the map f : X dom (/) 

X be continuous. 

Let Z C M", xq e Z, Z C dom (/). We define 

W^izo, f) = {z\ V„>o/"(z) e Z, hm /"(z) = zq}, 

n— >ctt 

^z(^Oi /) = {^\^ {^7i\ C Z a full backward orbit through z, such that 
lim Xn ^ zo}, 

71— '■ — OG 

W'{zoJ) - {z\ lim /"(z) = zo}, 

n — *oo 

W^{zo, f) — {z \ 3 {xn} a full backward orbit through z, such that 
lim Xn = zo}. 



n — ^ — 00 



// / is known from the context, then we will usually drop it and use W''^{zo), 
W^{z()) etc., instead. 



2.2 Projectivization of the dynamics 

The notation and setting is the one used in [AMj . Let f : X ^ X he a, 
diffeomorphism of a manifold X. 

Definition 3 Let zo E X . We say that zq is a hyperbolic fixed point for / iff 
f{zo) — Zq and Sp{Df{zQ)) D —9, where Df(zo) is the derivative of f at zq, 
Sp{A) denotes the spectrum of a square matrix A. 
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We denote the tangent bundle of X by TX and the differential of / by Df. 
From the dynamical system / : X ^ X we can derive a new dynamical system 
Pf : PX — > PX which is defined as follows. The space PX is the projective 
bundle associated to the tangent bundle of X, that is, the fiber bundle on X 
whose fiber over x € X is the projective space of T^X. That is, 

PX = P^X := {one-dimensional subspace of T^X}. 

xex xex 

For a submanifold S' in X by PS we will denote its projectivization, which 
is given by 

PS = {{x, [v]) ePX\xeS,veT^S\ {0}} 

It is easy to see that P{S) is a manifold and dim(P(5)) = 2dim(5') — 1. 

Define Pf to be the map induced from Df on PX, namely, Pf{x, [v\) := 
(/(a;), [Df{x) ■ v\) where ^ v G T^X, [v] is the subspace in T^X spanned by v 
and Df{x) : T^X Tff^-jX is the derivative of / at x. 

Let us identify X with the zero section of TX. We define the map tt : 
TX\X ^PXhyTT{{x,v)) = {x,[v])fovxGX,vGT^X. 

Let p e X be a hyperbolic fixed point of / and let TpX = Ep Q) Ep be the 
corresponding splitting of the tangent space into stable and unstable subspaces 
for Dfip). Define E; := 7r{E; \ {0}) and E^ := 7r(^« \ {0}). 

2.3 Dimension two 

Theorem 1 Let / : ^ be a diffeomorphism and let p be a hyperbolic 
fixed point with one dimensional stable and unstable manifolds. 
Then 

1. (p, Ep) is a hyperbolic fixed point for Pf such that 

. W^{ij>,E;),Pf)=PiW"ip,f)), dimW"{{p,E^),Pf) = 1, 

• W'{{p,E;),Pf) = {{z, M) e PX\z e W'{p,f), {z, [v]) i P{w\p,f))}, 

dimW%{p,E;),Pf) = 2. 

2. (p, Ep) is a hyperbolic fixed point for Pf such that 

. W'{{p,E;),Pf) = P{W'{p,f)), dimW'{{p,E;),Pf) = 1, 

• W"{{p, E;),Pf) = {{z, H) e PX\z e w-{p, /), {z, [v]) ^ P(W"(p, /))}, 

dimW-{{p,E;),Pf) = 2. 

Proof: It is easy to see that the stable and unstable sets are as stated. There 
remains for us to show the hyperbolicity, only. Consider first {p,Ep). Let us 
change the coordinate system in such that p = 0, Df{p) is diagonal 
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where |A| > 1 and < 1. Prom now on all considerations will be done in this 
coordinate frame. 

Around (p, £'"(p)) we will use the coordinate system (p iM? xM. ^ P(]R^) 
given by 

ip{z,v) = {z, [{l,v)]). 

It is easy to see that 

ip~^{z, [{V1,V2)]) = {z,V2/vi). 

In these coordinates {p^E^) is given by (0,0). Now we compute the lin- 
earization of Pf around (0, 0). 
We have 

Pf{z, v) = {f{z), [Df{z) ■ (1, vf]) = {Df{0) ■ z + o{\z\), [Df{z) ■ (1, vf\). 
Observe that 

(i\ A o\ (1 



jj, I \v 



\ + Oi\{z,v)\) 



.1^(0)^ + ^(0)2/ ^(0)^+1^(0)2/.' 



I^v+^{0)x+§^i0)y + o{\{z,v)\) 



Now we have to divide the second component of the above vector by the first 
one. Therefore 

V2/v,= (^^^v + ^{0)x+^iO)y + o{\{z,v)\)^ ■ X-'{1 + 0{\{z,v)\) = 

We have proved that the linearization of Pf at p = (0, 0) has the following form 

/A 

DPf{p,El)=[ M 

Vi0(o) iS(o) r 

Hence we see that {p, ) is a hyperbolic fixed point for Pf with one-dimensional 
unstable and two-dimensional stable manifolds. 

Let us consider now {p,Ep). This time we will use the coordinate system 
: M2 X R ^ P(IR2) given by 

= {z, [{v,l)]). 
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It is easy to see that 

(p~^{z,[{vi,V2)]) = {Z,vi/V2). 

Similar computations lead to the following formula for the linearization of Pf 
at {p,E;) 

/A 0\ 

DPf{p,E;)^i . 

V fi dxdy ^ ^ fi dy ^ ^ / 

I 

From Theorem [T] we obtain the following 

Remark 2 If f is as in Theorem[l\and W^lp, f) and W"{p, f) have a nonempty 
intersection, then 

• if W^lp, f) and W"{p,f) are tangent, then we obtain a heteroclinic con- 
nection from {p, Ep) to {p,Ep), 

• /) and W^{p,f) intersect transver sally, then we obtain a homo- 
clinic connection from (p, E^) to {p,Ep). 

3 Generic unfolding of quadratic tangency as 
the transversality question 

We assume that we have two curves in depending on some parameter a and 
given by u{a,t) = [fi{a,t), f2{a,t)) and s{a,t) = {gi{a,t),g2{a,t)). We are in- 
terested in establishing conditions, which will imply the existence of the generic 
unfolding of the quadratic tangency between them. Our goal is to formulate 
such conditions as the transversality question. This is Theorem 2.1 from |AM] . 
where it was stated without proof. 

Definition 4 JF^ Sec. 3.1] Let I,J,ZciR be intervals. Let u^,: I ^ and 
s^i : J for ji E Z be two smooth curves depending on ii in the smooth way, 

such that Ufi,g{tu) = Sf^„{ts) = qo and u and s are tangent at qq. 

Assume there exists a ^.-dependent coordinates in a neighborhood of q for p, 
close to fiQ, such that in these coordinates we can use xi (the first coordinate) 
as the parameter of our curves and the following holds 

s^ixi) = (a;i,0), 

Uf_,{xi) = {xi,axl-hb{^j.- Ha)) (2) 

where a ^ 0,b ^ 0. Then we say that the quadratic tangency of u and s unfolds 
generically. 

Remark 3 It is easy to see that in the above definition we can exchange the 
role of curves u and s using the coordinate transformation given by (j3fj,{xi, X2) = 
{xi,X2 - ax\ - b{^i - ^o))- 
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Remark 4 In the context of Def.^ = g{fi,xi) then instead of (0) it 

is enough to require 

5(/io,0) = 0, 1^(^0,0) = 0, 

OXi 

^W0)^0, 0(^0,0)^0. 
Proof: Observe that from our assumptions follows that 



Uf^{xi) = {xi,axl + b{fi - /io) + A{fj.,xi)), 



where 



a^Q, b^Q, 

A(mo,0) = 0, — (mo,0) = 0, —{^lo,0) = 0. 

Observe that through the parameter dependent coordinate chart (j)fj,{xi, X2) = 
(a:i, a;2 — A(^, xi)) we obtain (l2|). I 

Theorem 5 Let A C R be an interval. Let u : A X K ^ R2 and s : A X M ^ R 
he two -curves depending on parameter a G A. Let ao £ A be a parameter at 
which curves Ua and Sa are tangent. 

These curves have a quadratic tangency at qq which unfolds generically iff 
there exists tu,ts G M, such that surfaces 



E{u),E{s) : A X 



given by 



E(u){a,t) — ^,u{a,t), 
E{s){a,t) — I a,s{a,t), 



{a,t) 



du 

at 



ds , 



intersect transversally at the point E{u)[aQ,tu) — E{s)(ai^,ts). 
Before the proof we need one lemma. 

Lemma 6 Let s = (si, S2) : A x M ^ be a C'^-map, such that ^{oq, to) =/= 
for some (ao, to)- 

Then there exists a neighborhood A' of oq, an open set V , such that s{ao, to) € 
V, and a-dependent coordinates on V for a G A', i.e. (fia : V M.'^ for a G 
A', such that after a suitable reparameterization in these new coordinates the 
mapping s has locally the following form 



s(a,t) = (t,0). 
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Proof: Let us denote by {x,y) coordinates in R^. By taking a suitable param- 
eterization and shifting the coordinates origin and permuting, if necessary, the 
coordinates in we can assume that to = 0, s(ao,0) = and ^(ao,0) > 0. 
We can locally (in a suitable open set U) use x as the parameter of a curve 
s(a, •) for a S A C A. Therefore, we have 



s{a, t) = {t, S2{a, t)). 
Consider the following map ip : A x U ^AxM^ given by 
ip{a, X, y) = (a, x,y- §2(0, x)). 

We have 

D(f{a,x,y) = 



ds2 



1 



-^(a,x) 

oa ^ ^ ' ax 

Therefore, is a local diffeomorphism. Let open sets A', V be such that : 
A' X y ^ A' X is a diffeomorphism on the image and (ao, 0) G A' x V. The 
a-dependent coordinates on V are given by (f)a{x, y) = Lp{a, x, y). 

It is easy to see that in the new coordinates given by (pa the curve s has the 
following form 

s(a,t) = (t,0). 



Proof of Theorem [5} By Lemma [5] we can assume that the curve s is given 

by 

s{a,t) = (t,0). 

Observe that the transversality implies that in the neighborhood of the inter- 
section point the curve u can be represented as follows 

u(a,t) = {t,g2{a,t)), 

where 52 satisfies the following conditions 



<?2(ao,0) = 0, ^(ao,0)=0. 



(3) 



Now we will prove that the transversality of E(u) and E{s) is equivalent to 

^(ao,0)^0, 
oa 



d'92 



(ao,0) ^0. 



Observe that from Remark 2] it follows that the above conditions together with 
([3]) are equivalent to the generic unfolding of the quadratic tangency. 

Observe that ff(a,t) = (1,0) and ^{ao,t = 0) = (1,0) therefore in the 
neig hborhood of [ff(a,t)] we can use the second coordinate as a chart map in 
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In these coordinates we have 

E{s){a,t) = (a,i,0,0)^ 



E{u){a,t) 



(a, 32(0,*), ^{a,t)Y' . 



We have 



^ E{s){ao,t=Q) = span 



/I 







1 








^0 


0/ 



and 



-B(«)(ao,t=o) = span 



1 
1 

|f(ao,0) ^(ao,0) 

V|tK,0) ^(ao,0)/ 

From ([3]) it follows that the transversality question is equivalent to the following 
determinant being nonzero 



det 



/I 
1 



|^(ao,0) 



^KO)— (ao,0). 



\0 |t(ao,0) ^(acO)/ 
This finishes the proof. 



4 How to prove the homoclinic tangency us- 
ing the covering relations and the cone condi- 
tions? 

We assume that the reader is familiar with the following notions: h-sets, covering 
relations, cone conditions and horizontal and vertical disks as defined in [KWZ|, 

We consider a planar map fa'-M.'^D dom (/) ^ depending on the pa- 
rameter a, which has a hyperbolic fixed point pa with one-dimensional unstable 
and stable manifolds. Hence according to the setting from Section [2] we will 
work in four-dimensional space, using coordinates (a, x, y, v), where a is the pa- 
rameter, {x,y) e and v represents points in P(r(^ j,)R^), which we will call 
the tangential coordinate. We have the map (we abuse the notation for Pf) 

Pf{a, X, y, v) = (a, Pf{x, y, v)). 

In our method we need the following ingredients 
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the chain of covering relations 

No^N,^...^Nk, (4) 
such that the cone conditions are satisfied, 

{a,W'^{{pa, EpJ, Pfa)) as the horizontal disk in Nq satisfying the cone 
conditions, 

(a, W^{{pa, EpJ,Pfa)) as the vertical disk in Nk satisfying the cone con- 
ditions. 



If the above conditions are satisfied then from [Z] Theorem 7] it follows 
that (a, W^{{pa, Ep^), Pfa)) contains a horizontal disk satisfying the cone condi- 
tions in Nk. Hence (a, VK"((p„ £;^J, Pfa)) and (a, W'^Pa, E^J, Pfa)) intersect 
transversally in Nk, which by Theorem [5] implies that sets i¥'"((pa,Ep^),Pfa) 
and W'^{{pa, Ep^), Pfa) have a quadratic tangency which unfolds generically. 

Since the parameter a is not changing under Pf, apparently there is a prob- 
lem with the realization of the above scenario, because the covering relations 
together with the cone conditions imply the hyperbolicity. The essential point 
is that this remark is valid for closed loops of covering relations, but in our 
setting we just want a chain of covering relations going from one part of our 
phase space to another one. While constructing such chain we will arbitrarily 
decide whether we treat the parameter as the " unstable" /" stable" direction, by 
simply adjusting the sizes of h-sets in a-direction. 

To make our scheme to work we need: 

• to set the dimensions u, s in our h-sets Ni to be equal to 2, because this 
is the dimension of (a, W^{{pa, E^J,Pfa) and (a, W'{{pa, E^J, Pfa), 

• (a, W^iipa, Ep^), Pfa) should be a horizonal disk in Ng, hence we should 
treat a as one of unstable directions, 

• {a,W^ {{pa, Ep^), Pfa) should be a vertical disk in Nk, hence we should 
treat a as one of stable directions. 



4.1 The Lipschitz dependence of stable and unstable man- 
ifolds on parameters 

We need prove that (a, W'^{{p, Ep), Pfa) is a horizontal disk in a suitable h-set 
No, this requires at least the Lipschitz dependence on a of the invariant manifold. 
In [Z', Sec. 8.2] the Lipschitz dependence of (un)stable manifolds with respect 
to parameters with explicit and computable constants was discussed with the 
eye toward the computer assisted proofs. Here we will just recall (and refine a 
bit some estimates) these results in a form of the cone conditions. 

We will be using the norms for quadratic forms (identified in the sequel with 
symmetric matrices) which are defined by 

\B{u,v)\ < 
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For Euclidian norm we have 

||_B|| — max{|s| | s in an eigenvalue of _B}. 

Theorem 7 Assume that {N, Q) is an h-set in M"+* with cones and f\ : M"+* 
]^ti+s ^j^/j A e C, where C is a compact interval in the parameter space and Q 
has the form Q(x, y) = a{x) - I3(y) = Ya=i "'i^i ~ Y,i=i o-i+uuf ■ 

fx 

1. Assume that for the covering relation N N the cone condition is satisfied 
for all X e C. 



2. Let e > and A > be such that for all X G C and zi, Z2 € N holds 

Qihizi) - fx{z2)) - (1 + e)Q{zi - Z2) > A{zi - 22)'. 

3. Let 

d^zifx , 



M = max 
xec.zeN 



dz 



dX 



L 



\3\\ ■ max 



dX 



4. Let r > &e such that 



5. We define 



A - 2MT - LV^ > 0. 



Hall 



(5) 

(6) 
(7) 

(8) 
(9) 



Then the set Wf^{p\^ f\) for A G C can he parameterized as a vertical disk 
in C X N for the quadratic form Q{X, z) = 5Q{z) — X^ . 

Before the proof let us make two observations concerning constants A, e, T . 

Remark 8 The existence of A and e in {^i is a consequence of the cone condi- 
tion for covering relation N N . We would like to have as big A as possible. 
This forces e — > 0, but e is not used in the sequel. 

Remark 9 Since A > 0, therefore T in (0j always exists, but it is desirable to 
look for largest T possible. 

Proof of Theorem[7) We would like to obtain that for |Ai — A2I < — Z2\\ 
holds 

QifxAzi) - fxAz2)) > (1 + e)Q{z, - Z2). 
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Let i? be a unique symmetric form, such that B(u, u) — Q{u). Observe that 

Q(/ai (^i) - hjz2)) - (1 + e)Qiz, - Z2) = 
Q(/ai [zi) - fxAz2)) - (1 + e)Q(2i - Z2) + 
2B{fxAzi) - fxAz2)JxAz2) - fxAz2)) + Q{fxAz2) - fxAz2)). 

The first term in the above expression will be estimated using ([5]). 
For the third term we obtain 

Q{fxAz2) - fxAz2)) > ~P{7Ty{fxAz2) ~ fxAz2))) > 

2 

•(Ai-A2)2> 

2 

■T^\\zi - Z2f = ~LT^\\zi - Z2f. 

Finally, for the second term we have 

\B{fxAzi) - fxAz2)JxAz2) - .fxAz2))\ < 

•r||zi-z2||'. 

From the above computations and ([SHT]) we obtain the following 

QifxAzi) - fxAz2)) - (1 + e)Q{zi ~ z^) >{A- 2MT - LT^) ||zi - Z2II'. 

For r and (5 as in ([8]) and ([9]) it is proved in Lemma 23] that, for 
Ai, A2 e C, Ai 7^ A2 and Zi £ W^{pXi, fxj for i = 1, 2 holds 

SQ{zi ~ Z2) - (Ai - \2f < 0. 

Hence W^{px, fx) for A e C is vertical disk in C x TV for the quadratic form 
QiX, z) = SQiz) - X\ I 

Comments: 

• Since y, A) ~ Sa{x) — 5(3{y) — A^, hence due to the negative sign in 
front of A^ it follows that when trying to represent the stable manifold 
as the vertical disk in C x iV, we must treat the parameter as a 'stable' 
direction in an h-set. 

• For the unstable manifold we have to take the inverse map and we obtain 
different 5. Since taking the inverse involves changing the sign of Q we 
end up with the following quadratic form 

Q{x,y,\) = 5a{x) + \^ - 5(5{y). 

Looking at sign in front of A^ we see that A appears as an 'unstable' 
direction in C x A^. 



max 

xec 



dX 



(Z2) 



-\\B\\ max 

AeC,zGAf 



dT^yfx 



dX 



{z) 



max 



di^zjx 



(z) 



dTTzJx 



dX 



(z) 
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5 A toy example 



In this section we show how to construct the chain of covering relations ([4]) 
discussed in the first part of Section [4] for a special model map with a quadratic 
tangency which unfolds generically. Our intention is that the reader may see 
and appreciate some details of the method, which are later hardly visible when 
we report on the computer assisted proofs for the Henon map and the forced 
pendulum in Sections [H] and [71 respectively. 

We define a map / : ^ depending on the parameter a as follows: 

• (0, 0) is a hyperbolic fixed point and in a neighborhood of (0, 0) the map 
fa is linear 



where |A| > 1 and < 1, 

• in a neighborhood of the point (1,0) we have the homoclinic tangency for 
a = 0. We assume that fa acts as follows 



When compared with the full problem of establishing the existence of generic 
unfolding of homoclinic tangency for the Henon map and the forced damped 
pendulum, the above model map avoids the problems related to providing the 
explicit estimates on the dependence of on the parameter, which are given 
in Theorem [71 

5.1 Chain of covering relations 

Let p = (0,0), then = [(1,0)] and = [(0, 1)]. Our goal is to construct a 
chain of covering relations 'linking' {p,Ep) with (p,Ep). 
The beginning of the chain: 

Let us first see how the covering relations look in a neighborhood of {p, E^). 
We use the chart (a;, y, v, a), where v i-^ [(1, v)]. Therefore the map F — {Pfa, a) 
works as follows 



We define an h-set, so that {x, a) are the 'unstable' directions, {y, v) are 
the 'stable' ones. It is easy to see, that there exists a sequence of h-sets Ni = 
(ci + [-Xi,Xi]) X [-yi,yi] x [-Vi,v.i\ x [-ai.ai], with cq = 0, Cfc = 1 and such 
that 



fa{x,y) = (Aa;,^y) 



/a(l + x,y) ^ {x^ + y + a,l - x). 



F{x,y,v,a) = {\x,^y, {fi/\)v,a). 




Observe that the necessary conditions are 



Ci+i + [-Xi+i,Xi+i] C Xci + [-\X\xi,\X\x, 
\fi\yi < y,+i, 
\^l/X\vi < Vi+i, 
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The end of the chain: 

In a neighborhood of {p,Ep) the 'unstable' subspace is (a;, w), where w cor- 
responds to the [{w, 1)]. The map F works as follows 

F{x,y,w,a) = {\x,fiy, {\/^)w,a). 

It is easy to see, that there exists a sequence of /i-sets Mj = + [—Xi,Xi]) x 
[-yuVi] X [-Wi-,Wi] X [-ai,ai], such that 

co=0, c, = 1, M, ^M,_i^...^Mo. 

Observe that the necessary conditions are 

(Cj+i + C + [-\ii\xi,\ii\xi\, 

Vi+i > 

Wi+i < \X/lJ,\wi, 

Switching from unstable to the stable manifold. We want the covering 

relation Nk Mg, where Nk and Ms are as above. 

In Nk the nominally unstable directions are {x, a), while in Mg parameterized 
by {x,y,w,a), where w [{w, 1)] the 'unstable' direction is {x,w). We have 

F{1 + X, y, V, a) = (x^ + y + a,l — x, —{2x + v), a). 

We would like to have the following homotopy for covering relation 

Gt{l + x,y, V, a) = {a + t{x^ + y), 1 — tx, — 2x — tv, to). 

Therefore if wc want a covering relation Nk =^ Mg, whore Nk = (1, 0, 0, 0) + 
[-Xk,Xk] X [-yk.yk] X [-Vk,vk] X [-ak,ak] and = {0,1,0,0) + [-Xs,Xs] x 
[—ys,ys] X [—Ws,Wk] X [—ak,a.k], then we need to satisfy the following set of 
inequalities 

ak-xl-yk> Xs, 
2a;fe -Vk> Ws, 

Vs > Xk, 

as > ak- 

It is easy to see that this set of inequalities has a solution. For example if 
a/j = A < 1, then we can choose (for some small e > 0) 

Xk = A/2, 
yk = Xs = A/3, 
Vk = Ws = {l- e)A/2, 
ys = (0.5 + e)A, 
as = {l + e)A. 
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Observe that the expanding direction x is stretched in w-direction and a direc- 
tion is stretched across a;-direction (in the target set). 

It is clear that we can easily build the desired chain of covering relation of 
the form 

F F F F F F F 

5.2 Cone conditions 

It turns out that the cone conditions at the beginning and the end of the chain 
Q are relatively easy to satisfy in the situation when dynamics is linear (there 
is no parameter dependence) in a neighborhood of the hyperbolic fixed point, 
otherwise the issue becomes delicate see Theorem [7] in Section 14.11 

As a rule in this subsection the coordinate order for h-sets iV, Af will such 
that always the nominally unstable coordinates are written first. 

At the beginning of the chain For covering relations Ni Ni+i the 
expanding directions are (x,a). Let Ni be an h-set with Qi-cones given by 

Qat, (x, a, y , w) = a^a;^ + /3ia^ - 7i j/^ - (5iW^ . 

F 

The map F is linear for Ni iVi+i therefore it is enough to check whether 
QN,+AF{x^a,y,v)) > QNi{x,a,y,v). 

We have 

QNi+i {F{x, a,y,v))- Qn. {x,a,y,v) ^ 
ai+iX^x"^ + Pi+ia^ - 7i+i^^i/^ - {^i/XfSi+iv'^ - 
[aix^ + /Jjfl^ - 7iy^ - 5iv'^) = 
(ai+iA^ - ai)x'^ + (/3i+i - f3i)a^ + (7^ - ^^7j+i)?/^ + {5i ~ {^i/XfSi+i) 

Hence we will have cone conditions satisfied, when for example a^+i = ai, 
7i+i — 7i, 5i+i — 6i. The only strict requirement is 

A+i > A- 

At the end of the chain. It goes through as easily as the previous case. 
For covering relations Mi Mi-i the unstable directions are {x,w). We set 

QMi{x, w, y, a) = AiX^ + Biw'^ - dy^ - DiC? . 

As above it is enough to whether 

QMi^AF{x,w,y,a)) > QM,ix,w,y,a). 

We have 

QA'U-t{F{x,w,y,a)) - QMi{x,w,y,a) = 

(A,_iA2 - A,)x' + (B,_i(A//i)2 - lB,)w" + 

(a_i - aAi')y' + (A-i - A)a^ 
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We see that we need to have 

A-i > A. 
For the remaining coefficients wc can set 

= = Bi, Ci-i = Ci- 

Switching from the unstable to the stable manifold. Consider the 

p 

covering relation A''^ =^ Mg ■ This time we are in the nonhnear regime. We use 
on Nk coordinates (x, a, y, v) and on Ms coordinates [x, w, y, a) 

This means that in these coordinates (we remove the shifts of the origin of 
coordinate frame to (1, 0, 0, 0) and (0, 0, 1, 0) for and Mg respectively) 



F{x, a, y, v) = {x'^ + y + a, —2x — v, —x, a). 



We set 



We have 



QN{x,a,y,v) 
QM{x,w,y,a) 



= ax^ + Pa^ — — Sv"^, 
= Ax^ +Bw^ -Cy'^ -Da^. 



DF{x,a,y,v) 



2a; 1 1 
-2 0-1 

-10 

10 



After some computations we obtain 

'AAx^ + {4B-C-a) 

Q = [DFfQMDF -Qn= 

2B 



2Ax 
A-D-13 
A 




2Ax 
A 

^ + 7 




2B 




B + 6 



We are interested whether Q is positive definite. Observe that since the positive 
definiteness is an open condition we will show that Q is positive definite for a; = 
and then we will know that the same holds for \x\ small. 
After we set .t = we obtain 



Q = 



It is easy to see that after rearrangement of coordinates the question is reduced 
to the positive definiteness of the following two matrices 



-C-a 








2B 





A-D-p 


A 








A 


A + 7 





2B 








B + 5 



Qi 



AB-C 
2B 



a 2B 
B + 5 



Q2 



A-D-(3 A 
A A + 7 
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For example, we can set A = B = C = 1 and D = \, then we obtain 

Q2 = 



3-a 2 
2 1 + 5 



1/2-/3 1 
1 1 + 7 



It is now easy to see that we get what we want, when we set 

a = 1, (5 > 1, 
/3 = l/4, 7>3. 

Summarizing we obtained 

QN{x,a,y,v) = + / A - Ay^ - 2v'^ , 
QM{x,w,y,a) = + - y'^ ~ /2. 



6 Application to the Henon map. 

In this section we will show that the method introduced in the previous sections 
can be successfully applied to a specific system. Let us consider the Henon map 

Ha,b{x,y) ^ (a- x'^ + by,x). (10) 

The following theorems have been proven in [AMj . 

Theorem 10 fAM\. Thm.1.1] Th ere exists an open neighborhood B of param- 
eter value b — 0.3 such that for each parameter b d B there exists a parameter 
value 

a e [1.392419807915,1.392419807931] 

such that the Henon map has a quadratic tangency unfolding generically for the 
fixed point 

Xa,b = ya,b = (b- y/{b-l)^ +4:0-1 



Theorem 11 fJM Thm.1.2] There exists an open neighborhood B of parame- 
ter value b = —0.3 such that for each parameter b d B there exists a parameter 
value 

a e [1.314527109319, 1.314527109334] (11) 

such that the Henon map has a quadratic tangency unfolding generically for the 
fixed point 



1 



Xa,b = yaM ^ ^{b-iy^ + Aa - 1 

The numerical evidence of the existence of a homoclinic tangency for some 
parameter values (a, b) « (1.3145271093265, —0.3) is shown on Fig. [T]- see also 
[Ml Fig. 1.1] 

The main motivation for us to study the existence of homoclinic tangencies 
for the Henon map was to verify if our method could work in this relatively easy 
example. 



17 



-2-1 1 2 



Figure 1: Parts of the unstable and stable manifolds of the Henon map at the 
fixed point for b = -0.3 and a = ao = 1.3145271093265. 

Denote by 

ao = 1.3145271093265 

a center of the interval pT|) . The aim of this section is to prove the following 
theorem 

Theorem 12 There exists an open neighborhood B of the parameter value b = 
—0.3 such that for each b G B there is a parameter a e ao + [—10^^, 10"^] such 
that the Henon map (|10p Ha.t has a quadratic homoclinic tangency unfolding 
generically for the fixed point 

Xa,b = Va.b ^\{^- +4a - l) w -1.9679632427827796. (12) 

The authors in |AMj report relatively long computational time (approxi- 
mately 100 minutes) for the computer assisted proof of [AMI Thm.1.2]. It turns 
out that using our method the verification of necessary inequalities in the com- 
puter assisted proof of Theorem [T2l has been completed in 0.2 second on the 
Intel Xeon 5160, 3GHz processor. This very good efficiency allowed us to apply 
the method to a map coming from an ODE - a suitable Poincare map in the 
forced-damped pendulum. The details will be given in the next section. 

Computer assisted proof of Theorem [T^] will be presented in the following 
subsections in which we verify: 

• the existence of a heteroclinic chain of covering relations for PH , 

• the cone conditions along the above mentioned heteroclinic chain, 
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• the cone conditions at the beginning and at the end of the chain which 
allow us to parameterize the center-unstable and center-stable manifolds 
as a horizontal or vertical discs, respectively, in proper h-sets. 



6.1 The existence of a heteroclinic chain of covering rela- 
tions for PH. 

In Theorem [TH we chose the center of the interval (HH) oo = 1.3145271093265 
as a good candidate for the homoclinic tangency parameter corresponding to 
60 = -0.3. 

In order to define the sets which will appear in the heteroclinic chain, we 
need to set a local chart for the manifold PR^. It turns out, that it is enough 
for our purpose to use the parametrization 

: X (0, tt) 3 {x, y, t) ^ (x, y, [(cos(i), sin(t))]) G PR^ (13) 

This parametrization excludes one point on the manifold, but trough our com- 
putations this point does not appear either as an argument or a value of PHa- 
All the sets will be expressed using coordinates {x, y, t, a). 
Put 

PH{x,y,t,a) - {{^-^oPHa,bo°i'){x,y,t),a) (14) 
In the sequel by tt^ we will denote a projection onto tangent coordinate, i.e. 

Trt{x,y, [u]) ^ 

\m\ 

With some abuse of notation we will use the same symbols for the projections 

u 

TTt{x,y,[u],a) = —— or 7rt(a;, y, t, a) = (cos(i), sin(i)) 

\m\ 

but it is clear from the list of arguments which projection has to be used. In 
each case the value of tti is a vector. 

To simplify the notation we will use zq — {xa,yo) = (a^ao.feoi UaoM) defined 
in p^ . Since we always will have the fixed value of the parameter 60 = —0.3 
we will write Ha instead of i?a,6o the sequel. Let uq and sq be normalized 
with respect to the Euclidean norm eigenvectors of DHao{zo) given explicitly 

by 



uo = (Z£±v^^£) ^ (0.9680131177714217873,0.250899589123719882), 

\\{-xo + \/x:l+boA)\\ 

So = (-^»-Vfp^° ii) ^ -(0.07752307795993337433,0.996990557820693689) 

(15) 

and let M = [uq , s'^] be a matrix of eigenvectors. Put 

zi=zo + 0.0001993152279412426wo + 2.50404 • 10""so- 
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The above point has been chosen as a good approximation of the homochnic 
tangency point for Hag ■ Namely, we have 



|ji?a"o'(^i)-^o|| <5.2•10-^ 
lK'(^l)-^o||<1.2•10-^ 
{iTt {PH^tizi, [uo]))) - (-4.71 • 10-^ 0.999999847). 



(16) 



We see that the unstable direction uq at zi is mapped under 14-th iterate of 
Hag very close to the stable direction sq at the point Hag{zi). Let us underline 
that to get the estimation we need to compute DH^"^ with at least long 
double precision of the floating point arithmetics. 

The points on the trajectory of zi will be the centers of the sets which appear 
in the computer assisted proof. Put 



Co = 

Cl = 

Ci+l = 

Cl5 = 

For further use we set also 

Zi+l = 

2^15 = 



(■(/' ^(zo,uo),ao), 
(V'"^(zi,uo),ao), 
PH\ci), for i = 1, 

(V'~^(2o,so),ao)- 



Zq. 



for 1 = 1, 



13, 



13, 



Some of these points are shown in Fig. [T] 

Now, we have to chose an approximate stable and unstable directions at Ci, 
i = 0, . . . , 15. On each set centered at Ci, the coordinate system will be given 
by a matrix 

"(m,)i (s,)i 0" 

iu^)2 is^)2 

10 
1 



M, 



(17) 



where ug, sq are given by (jlSp and Ui and Si are computed as follows. 

• Put Ul5 = Uo, Sl5 = So, Ul = Uo, si = So- 

• For i = 2, . . . , 8 we set 



Ui = TTtiCi), 

Si = Trt{PH~'^{Zi+l,TTt{Ci+l) 



-))• 



(18) 
(19) 



From numerical simulations we get, that between the points cs and cg the 
role of the tangent coordinate is changing from contracting to expanding. 
Therefore, the unstable direction propagates very well on these sets just 
by (fT5)) for i = 2, . . . , 8. For the inverse map, the unstable directions Ui 
become repelling. Therefore the preimage of an orthogonal direction to 
Ui+i becomes a good enough approximation for our method of the stable 
direction at Zi. 
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i 


10^ -{dih 


10^ • {d,h 


10^ • (d,)3 


lO*^ • {di)4 




unstable dir. 


stable dir. 


tangent dir. 


parameter 





7 


1 


2 


(1.01)« 


1 


1 


1 


2 


(1.01)^ 


2 


1 


1 


2 


(1.01)" 


3 


1 


1 


2 


(l.Olf 


4 


1 


1 


2 


(1.01)^ 


5 


1 


1 


2 


(l.Olf 


6 


1 


1 


2 


(1.01)^ 


7 


1 


1 


2 


1.01 


8 


1 


1 


2 


1 


9 


0.5 


1.25 


0.25 


1.01 


10 


0.75 


1.25 


0.25 


(1.01)^ 


11 


1 


1.25 


0.25 


(1.01)^ 


12 


1 


1.25 


0.25 


(1.01)4 


13 


1 


1.25 


0.25 


(1.01)" 


14 


1 


1.25 


0.25 


(1.01)« 


15 


1 


2 


0.25 


(1.01)^ 



Table 1: Diameters of the h-sets in the heteroclinic sequence for the Henon map. 
The diameters in the table are scaled by the factor 10^. On h-sets Nq, . . . , Ns 
the expanding directions are u and a coordinates, while on iVg,...,iVi5 the 
expanding directions are u and t coordinates. 

• According to a good choice of the homoclinic point ([TB|) . the tangent 
coordinate at Ci become a good approximation of the stable direction on 
sets centered at cg, . . . , C14. Therefore, we can set 

Si = 7rt(ci), for i = 9, . . . , 14, 

Ug = TTt{PHao{z8,[s8])), 

Ui+1 = T:t{PHaoizi,[ui])), for i 9, . . . , 13. 

There remains for us to set the sizes of the sets. We define the h-sets by 

N,=c,+M,-id,-[~l,l]), (20) 

where the diameters di are listed in Table [T] and chosen from numerical exper- 
iments as well working. On h-sets A'o, . . . , iVg the expanding directions are u 
and a coordinates, while on iVg, . . . , Ni^ the expanding directions are u and t 
coordinates. 

Let us comment briefly about the choice of the sizes of the sets presented in 
Table [T] For the sets Nq, . . . , Ns the tangent coordinate is chosen as a contract- 
ing direction. Therefore the set in this direction must be large enough (here 
2 • 10^^) to be able to easy verify that we have contraction on this variable. 
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On the other hand, for the sets Ng, . . . , Ni^ we have a strong expansion on this 
coordinate. Therefore the diameter of these sets on this coordinate is smaller 
here because it is easier to verify that the image of some walls of the previous 
set in the chain is outside of a small set than inside of it. 

Moreover, the tangent coordinate of Ng has to be covered by the unstable 
coordinate of iVg. This expansion is weak and forces decreasing of the size Ng 
on the tangent coordinate. 

The unstable and stable direction are comparable for almost all sets in the 
chain except Nq, Ng and A^is. For Nq we have large size in unstable direction. 
This is due to the fact that A^i is in some distance from zq. The size 7 • 10~^ is 
necessary to reach the set iVi from Nq . 

Similar reasoning apply to -/V15. The stable size must be large enough so 
that the image of N14 is captured in this direction. 

On the set Ng we have change of dynamics. The unstable size of A'g is 
smaller since on Ns the parameter plays a role of expanding direction which 
covers unstable direction on Ng. Since this expansion is weak, we must decrease 
the size of A^g . On the other hand the image of Ng vary with the parameter and 
it is larger in the stable direction of A^g. This forces to enlarge a little bit the 
set A^g in the stable direction. 

Now we can state the first numerical lemma. 

Lemma 13 The following covering relations hold true 



Proof: The assertion of the above lemma has been verified using the inter- 
val arithmetics [M] and the algorithms for verifying the existence of covering 
relations described in [WZ| . We were able to verify the necessary inequalities 
on each wall of A^,, i = 0, . . . , 15 without any subdivision of the sets, so the 



6.2 The cone conditions along the heteroclinic chain 

In this section we will show that the cone conditions are satisfied for the sequence 
of covering relations from Lemma 1131 We have two results which give us a 
numerical method for verifying the cone conditions. 

Lemma 14 \KWZ[ Lemma 6] Let {N,Qn) and {M,Qm) be h-sets with cones 

in R" and let f : N M" be such that N J=> M. Let [Df{N)]i denote the 
interval enclosure of the set of matrices Df[N). Lf the interval matrix 



is positive definite then the cone conditions are satisfied for the covering relation 




PH 



PH 



PH 



computational time was less than one second. 



V = [Df{N)f,QM[Df{N)]i - Qm 



N 




M. 
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Let A = Ac + [—1,1] ■ Aq he an interval matrix, where Ac,Ao G R"-^" are 
real and symmetric. For z € M" by A(z) we will denote a diagonal matrix with 
Zi's at the diagonal. 

Lemma 15 |^ Theorem 2] The interval matrix A = Ac + [— 1,1]- Aq is positive 
definite if and only if for each sequence z G { — 1, 1}" the matrix 



is positive definite. 

In the light of the above lemma to verify that a symmetric interval matrix 
is positive definite it is enough to verify if 2"^^ real matrices A^ are positive 
definite. 

In order to verify the cone conditions we have to define quadratic forms on 
the sets Ni, z = 0, . . . , 15. Denote by 



an approximate eigenvalues of DHa^ at zq. Recall, by A(pi, . . . ,p„) we denote 
a diagonal matrix with piS at the diagonal. For i — 0, . . . , 15 we define the 
quadratic form on the h-set Ni by 



where the coefficients are listed in Tabled] The quadratic forms Qi are defined 

in the coordinate systems given by matrixes Mi pT]) used to define the h-sets Ni, 

i = 0, . . . , 15. These matrices have normalized columns. In these coordinates 

the sets are given by Ni = [{di)i, . . . , {di)4] ■ [—1, 1], i = 0, . . . , 15 - see (pOj) 

and Table [TJ Let us comment briefly about the choice of these coefficients. 

Assume that we would like to verify the cone conditions for the covering relation 
f 

N M. Assume that / is linear and in some coordinate systems on N and M 
it is given by / = A(Ai, . . . , Afc). In general case, we usually have Df{N) close 
to a diagonal matrix, but the arguments apply. Assume also that the quadratic 
forms on both sets N and M are diagonal and given by 



According to Lemma [TU the cone conditions will be satisfied if the interval 
matrix 



V = [Df{N)]] ■ Qm ■ [Df{N)]j - = A(A2< - af , . . . , Afaf - a^) 



is positive definite. We see, that if ' | = 1 and some |Ai| ^ 1 then the 
corresponding coefficient in the matrix V becomes very large while for |Aj| <^ 1 
the corresponding coefficient in V is close to 1. In general (nonlinear) case 
we have some nonzero intervals off the diagonal of V. Therefore, from the 



A^^Ac- A{z)AoA{z) 



A = 3.858169402, ^ = 0.07775708341 



(21) 



Qi= A {{pi)l, {pi)2, {Pi)3, {Pi)4) 



Qn 
Qm 
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i 


(k)i 
unstable dir. 


{Pi)2 

stable dir. 


{Pi)3 

tangent dir. 


parameter 





3/A" 


-/^^ 


-WA)" 


2(1.5)-» 


1 


1/A" 


-0.1 


-0.5 


2(1.5)-^ 


2 


1/A" 


-0.1 


-1 


2(1.5)-" 


3 


1/A" 


-0.1 


-1 


2(1.5)-" 


4 


1/A" 


-0.1 


-1 


2(1.5)-^ 


5 


1/A" 


-0.1 


-1 


2(1.5)-^ 


6 


1/A" 


-0.1 


-1 


2(1.5)-" 


7 


1/A" 


-0.1 


-1 


2(1.5)-^ 


8 


G.5/A" 


-1 


-1 


2 


9 


100/A" 


-0.1 


100(Af/A)" 


-2 


10 


40/ A" 


-0.1 




-2(1.5)-^ 


11 


10/A" 


-0.1 




-2(1.5)-" 


12 


1/A" 


-0.1 




-2(1.5)-^ 


13 


1/A^ 


-0.1 


(a^/a)^ 


-2(1.5)-^ 


14 


1/A" 


-0.1 


WA)^ 


-2(1.5)-" 


15 


G.3/A" 


-0.1 


(a^/a)" 


-2(1.5)-" 



Table 2: Coef&cients of the quadratic forms on the sets Nq, . . . ,Ni^, where 
A = 3.858169402, = 0.07775708341 are approximate eigenvalues of DHao{zo). 



computational point of view, it is better to make the matrix V somehow uniform, 
i.e. such that the coefficients on the diagonal are of the same magnitude. This 
can be achieved by setting coefficients Ui « A"^ for these i such that \Xi \ 3> 1. 

This can be seen in the first column of the Table [2] and the second part of 
the third column. 

The coefficients for iVg are chosen to be able to switch from the unstable to 
the stable manifolds. 

Notice also a different choice of the coefficients for the quadratic form Qo 
in Nq. The reason is that we have to prove that the center-unstable manifold 
at the fixed point cq is a horizontal disk in A'o satisfying the cone conditions. 
Therefore we need to apply the method described in Theorem[7]to the inverse of 
PH. In this case, the coefficients in the second and the third column correspond 
to expanding directions with approximate eigenvalues l//i and A//i, both greater 
than 1. 

We have 

Lemma 16 The cone conditions are satisfied for the sequence of covering rela- 
tions 

Proof: Observe that verification of the cone conditions for these covering 
relations require computation of the DPH which involves the second order 
derivatives of H . In the computer assisted proof we computed an enclosure for 
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DPH{Ni), i = 0, 1 . . . , 14 using whole Ni without subdivision as an initial con- 
dition of the routine which computes DPH. Then we apphed Lemma 1 141 and 
Lemma [15] to prove that the cone conditions are satisfied. The C++ program 
which verifies the assertion executes within less than one second on a laptop- type 
computer. I 



6.3 Parameterization of center-unstable and center-stable 
manifolds at cq and C15, respectively. 

In this section we will use the method described in Theorem [7] in order to 
parameterize the center-unstable and center-stable manifolds as a horizontal and 
vertical discs satisfying the cone conditions in Nq and iVis, respectively. This, 
together with Lemma [13] and Lemma [TH] will give us a proof of Theorem [T^] 
Put 

No = (2o,M) + M- ((rfo)i,(do)2,(rfo)3), 

7Vi5 = (2o,N) + M-((di5)i,(di5)2,(dl5)3), 

where M is a 3 x 3 minor of the matrix A/15 = Mq (see ([T7|) after removing the 
last column and the last row. The coefficients {di)j are listed in the Table [1] 
Geometrically, these sets are just projection onto {x, y, t) coordinates of Nq and 
7Vi5, respectively. 

The set iVis is a three-dimensional h-set with two expanding directions (cor- 
responding to eigenvalues A and ^) and one nominally stable direction (cor- 
responding to the eigenvalue ^), where A,/i are eigenvalues of DHao{zQ) - see 

m- _ 

On the set A^o we will compute the inverse map of PHa so that the role of 
nominally stable and nominally unstable directions interchange. Hence, the set 
A^o has two nominally unstable directions (with eigenvalues and A//i) and 
one nominally stable direction (with eigenvalue A~^). 

On both sets we set quadratic forms defining the cones to be equal to 

Qo = -A((po)i,(Po)2,(po)3), (22) 

= A((pi5)i,(pi5)2,(pi5)3), (23) 

where the coefficients are listed in the Table [2] 

Lemma 17 The center-stable manifold of PH at C15 can be parameterized as a 
vertical disk in Ni^ satisfying the cone conditions with respect to the quadratic 
form Q15. 

Proof: We use Theorem [7] Let tTq denote a projection onto the parameter 
coordinate. With a computer assistance we verified that for a G iTaiNi^) holds 

PH 

Ni5 ^ iVi5 
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and the cone conditions are satisfied with the quadratic form Q15. Hence, the 
first assumption of Theorem [7| is fulfilled. 

Then we computed the constants A, M , L which appear in assumptions 2-3 
of Theorem [7] and we obtained 

A > 0.099394300936541294, 
M < 0.084042214456891598, 
L < 0.0070394636406844067. 

Hence constant F from assumption 4 of Theorem [7| can be chosen to be equal to 
r = 0.57737423322563175. With this T the coefficient S defined in ^ is equal 
to 

S^---- < 16.54078540195168. 

Il"ll 

where a appears in the decomposition of the quadratic form Qi^{x,y,t) = 
a{x,t) - (i\y) and ||a|| = max{0.3/A2, [y./\f] = O.S/A^ - see TableU 

From Theorem [7] it follows that the center-stable manifold of PH at C15 can 
be parameterized as a vertical disk in iVis satisfying the cone conditions with 
respect to the quadratic form 

Qib{x, y, t, a) = SQi5{x, y, z)~a^. 

Recall that by Q15 we denote the quadratic form on N15. To finish the proof 
let us observe that 

Qi5{x, y, t, a) = oT^ix, y, t) - 2{1.5)-''a^ 
see (I2S1) and Tabled] Moreover, we have 2(1. 5)"^^ > 1. This shows that 

SQisix, y, t,a) (^Qi5{x, y, t) - 2(1.5)"^a^^ < Qizix, y, t, a). 

Therefore the center-stable manifold of PH at C15 is a vertical disk in iVis 
satisfying the cone condition for the quadratic form Qis. I 

We have a similar lemma about the parameterization of the center-unstable 
manifold of PH at cq as a horizontal disk in TVq. 

Lemma 18 The center-unstable manifold of PH at cq can he parameterized 
as a horizontal disk in Nq satisfying the cone conditions with respect to the 
quadratic form Qo ■ 

Proof: We will proceed as in Lemma [T71 but for the map PH^^. With a 
computer assistance we verified that a e Tra{No) holds 

— PH^^ ~ 
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and the cone conditions are satisfied. We computed the constants which appear 
in dlHS]) and we got 



A 


> 


0.1877584261322994, 


M 


< 


0.2795983187542756, 


L 


< 


0.015049353557694945, 


r 




0.33278415598142302, 


5 


< 


18.316620936531205. 



Hence, the center-stable manifold for PH ^ at cq is a vertical disk in Nq satis- 
fying the cone condition for the quadratic form 



Qoix,y,t,a) = 5Qo{x,y,t) 



,2 



a 



From (EH) and Tabled] we have 

Qq{x, y, t, a) = -Qo{x, y, t, a) + 2(1.5)~*a^. 
Since 2(1.5)~*(5 > 1 we have 

-6Qo{x, y, t,a) = 5 {Qo{x, y, t, a) - 2(1.5)"^a^) < Qo(x, y, t, a) 

and the center-stable manifold of PH~^ at cq can be parameterized as a vertical 
disk in Nq satisfying the cone conditions with respect to the quadratic form — Qo- 
This means that the center-unstable manifold of PH can be parameterized as a 
horizontal disk in Nq satisfying the cone conditions with respect to the quadratic 
form Qq. I 

Proof of Theorem I12t We obtain our conclusion from Lemmas [HI [T7] and 
[T8l combined with the discussion of the strategy of the proof in the first part of 
Section m I 



7 Application to the forced damped pendulum 
equation. 

Let us consider an equation for the forced damped pendulum motion 

X + l3x + sm{x) = cos{t). (24) 

The equation l|24p is a well known example of an equation which exhibits chaotic 
dynamics see |GHj and references given there. In |BCGHj it has been proven 
(a computer assisted proof) that for the parameter /3 = 0.1 the 27r-time map 
is semiconjugated to the full shift on three symbols on some compact invariant 
set. 
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The equation ([M)) defines a flow on x S'-'^, where is a unit circle. Let 
us define the Poincare map : ^ by 

Tp{x,i) ^ix{2Tr),x{2Tr), (25) 

where x{t) is a solution of ([M)) with the parameter /?. 

The aim of this section is to prove the following theorem. 

Theorem 19 For all parameter values 

(3eB = 0.247133729485+ [-1, 1] • 1.2 • 10^1° 

there exists a hyperbolic fixed point for Tp corresponding to a 2tt periodic solution 
of {2Ji}jj . Moreover, there exists a parameter value [3 ^ B such that the map Tf} 
has a quadratic homoclinic tangency unfolding generically for that fixed point. 




2468 10 12 345678 



Figure 2: Left: parts of the unstable and stable manifolds of a fixed point of T^, 
for (3 — 0.247133729485. Right: Periodic orbit and a homoclinic orbit to the 
periodic orbit for pi)) projected onto {x, x) coordinates. This homoclinic orbit 
corresponds to the quadratic tangency of invariant manifolds. The left picture 
suggests also the existence of transversal homoclinic points. 

The proof of Theorem [T9l uses the same method as for the Henon map case. 
Since all the details has been discussed in Section [6] we will give here only the 
definition of the sets, quadratic forms and we will state the numerical lemmas. 

The main difference between the Henon map case and that of the forced 
damped pendulum equation is how we compute derivatives of a map up to the 
second order, required in the proof. Since the Henon map is given explicitly these 
derivatives can be computed by hand or using the Automatic Differentiation 
tools [G]. For the forced damped pendulum equation we used the C"'-Lohner 
algorithm presented in [WZ2j . This algorithm allows us to integrate efficiently 
the variational equations for ODE's and compute partial derivatives of Poincare 
maps. 



28 



7.1 Heteroclinic chain of covering relations. 

Let PTp : PK^ — > PR^ denote the map induced by on the projective bundle 
and let -0 be a local parameterization of the space PR^ as defined in . Let 

PT : PR2 X R 9 {x,x,t,/3) = {{t/r^ o PP^ o iP){x,i,t), (3) e PR^ x R. 

Let 

/3o = 0.247133729485, 

zq = (2.6410109874338904,0.063471204982120187), 

uo = (0.76818278871270265,0.64023058590290362), (26) 

So = (0.67655372773981033,-0.73639327365298942), (27) 

A = 211.83022271012155, (28) 

= 0.00099918347695168025. (29) 

From numerical experiments we got that zo (see Fig. [2]) is an approximate fixed 
point for TfSg with approximate eigenvalues A, fi and approximate eigenvectors 
Mo, So, respectively. This point has been found by the standard Newton method. 
Let us denote 

Co = V'"^((2o,wo),/3o), 
Zl = zo + 1.41442890240556- 10-%, ci = ^"^((^i, "o), /3o) 
Zi = P/3„(2i-i), for z = 2,3,4, 5, = PT(ci-i), for i = 2, 3, 4, 5, 

zs = Zo + 19.0992395815- 10-«so, cg = ^"^((^8, so), /3o) 
Zi = P3"^^(zi+i) for I = 7, 6, = PT-i(ci_i) for I = 7, 6, 

= Zq, Cg = V"^((20,So),/3o)- 

Some of these points are shown in Fig. [2l The points zi and zg are chosen close 
to a heteroclinic trajectory for Tp^. Numerical simulation shows that 

PT-\ce) - C5 « (3.42 - lO^^^ 6.34 - 10"l^ -2.32 - 10"^ 0). (30) 

In fact, we observe that between cs and cg the role of the tangent direction 
changes and at cq it becomes repelling. 

The reason for which we define two points zi and zg and compute their 
forward and backward trajectory, respectively, is due to the numerical problems 
with forward propagation of tangent coordinate at zg, zy, even in simulation 
only. Recall, the tangent coordinate t becomes strongly repelling at zg and 
Zl with eigenvalue of the order 10^. Therefore, it is easier to propagate the 
preimage of the stable direction at zg which is attracting at zj and zg for the 
inverse map. 

Approximate estimation ([50)1 shows us that the points c^, i ~ 1, ... ,8 are 
close to the possible existing heteroclinic trajectory for PT. 

Now we will define the coordinate systems of the h-sets centered at c^'s. Let 
Mi be the matrix of coordinate system of the set centered at Ci . We assume Mi 
has the form (fT7|l . These matrices are computed as follows 
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i 




lO^'^ • {di)2 


10« • {d,)3 


lO^'^ • (d,)4 




unstable dir. 


stable dir. 


tangent dir. 


parameter 





5 


0.2 


0.26 


1.2(1.01)^ 


1 


0.4 


6 


2 


1.2(1.01)" 


2 


0.4 


6 


4 


1.2(1.01)^ 


3 


0.6 


6 


2 


1.2(1.01)^ 


4 


5 


2 


1 


1.2(1.01) 


5 


27 


2 


1 


1.2 


6 


0.4 


15 


0.4 


1.2(1.01) 


7 


0.4 


8 


0.4 


1.2(1.01)^ 


8 


0.4 


10 


1.2 


1.2(1.01)^ 


9 


0.4 


5 


0.2 


1.2(1.01)" 



Table 3: Diameters of the h-sets in the heteroclinic sequence for the PT map. 
The h-sets Nq, . . . , TVs have two unstable directions given by Ui and the param- 
eter. The h-sets Nq, . . . , Ng have two unstable directions given by Ui and the 
tangent coordinate. 



• put ug = uq, Ui = uo, 1*8 = uq and sg = sq, Si = Sq, ss = So, where 
Mo, Sq are defined by 

• for z = 2, 3, 4, 5 we set 



Ui 
Si 
S5 



= 7rt(ci), 

= 7rt(PTg;^^(z,+i,7rt(ci+i)^)),i = 2,3,4, 



*5 : 



• for i = 6, 7 we set 



7I't(Ci), 



As in the case of the Henon map the h-sets Nq, . . . ,Ng for the map PT will 
be defined by formula ((20|) with matrices Mi as above and the diameters given 
in Table H 

The h-sets Nq, . . . , have two unstable directions given by Ui and the 
parameter. The h-sets N^, . . . , Ng have two unstable directions given by Ui and 
the tangent coordinate. 

Now, we can state the following numerical lemma 

Lemma 20 The map PT is well defined and continuous on IJ^^q Ni . Moreover, 
the following covering relations hold true 



PT 



No. 



(31) 
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Proof: In the computer assisted proof of the above lemma we used the C^- 
Lohner algorithm [Z] and the CAPD library [CAPDj in order to integrate the 
variational equations for and to compute the map PT. We used the Taylor 
method of the order 20 with a variable time step. To verify the inequalities 
required for the covering relations we subdivided the boundary of each h-set with 
a grid depending on the set under consideration. The total number of boxes 
we used is 5546. The C++ program which verifies the existence of covering 
relations (|3ip executes within 18 seconds on a computer with the Intel Xeon 
5160, SGHz processor. I 

7.2 The cone conditions along the heterochnic chain of 
covering relations for the map PT. 

Recall, that for a p G M" by A(p) we denoted a diagonal matrix with pi's on 
the diagonal. 

For i = 0, . . . , 9 we define the quadratic form on the h-set Ni by 



where the coefficients are listed in Table SI 

Let us comment about the choice of these coefficients. In the example pre- 
sented in Section [5] there is no dependency on the parameter, hence the cone 
conditions are easily achievable. In general case this dependency has huge in- 
fluence for the choice of the parameters in quadratic forms corresponding to the 
parameter variable. In fact, this sometimes forces large scaling (not only by 
some small factor), like in the Henon map case). 

The other constraints on the coefficients are related to the parameterization 
of center-unstable and center-stable manifolds in the first and last sets in the 
heteroclinic chain of covering relations. 

The main constraint, however, appears for the covering relation in the switch 
between manifolds. Here it is necessary to set relatively large coefficients cor- 
responding to both unstable variables in the set after this switch and small 
coefficient for the stable variable in the main phase space. This can be seen in 
Tables m and 0] for the sets A'g and A^g, respectively. In the next sets we can use 
the hyperobolicity of the map to make these coefficients more uniform and to 
reach constraints at the begin and at the end of the chain. 

For i = 2, . . . ,9 the coordinate systems on Ni are given by the matrices Mi 
used to define the sets Ni, respectively. 

We have the following numerical result 

Lemma 21 The cone conditions are satisfied for the sequence of covering rela- 
tions 



Proof: By Lemma [14] it is enough to verify that for i = 0, . . . , 8 the interval 
matrix 



Qi = A ((Pj)i, (Pj)2, (Pi)3, {Pz)4) 




V, = [DPT{N)]jQN,^,[DPTiN)]i - Qjv, 
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i 


iP^)l 




(k)3 






unstable dir. 


stable dir. 


tangent dir. 


parameter 





80/A" 






(1.1)-^ 


1 


1/A" 


-0.01 


-10-^ 


(1.1)-* 


2 




-1 


-10-'^ 


(1.1)-^ 


3 




-1 


-10-^ 


(1.1)- 


4 


1/A^ 


-1 


-10-^ 


(1.1)-^ 


5 


10/A" 


-1 


-10-'^ 


3 


6 


1000/A^ 


-10-4 




-(1.1)-^ 


7 


1/A^ 


-1 




-(1.1)- 


8 


1/A^ 


-1 


WA)^ 


-(1.1)-^ 


9 


1/A^ 


-1 


WA)^ 


-(1.1)-* 



Table 4: CoefBcients of the quadratic forms on the sets iVo, . . . , iVg, where A, fj, 
are approximate eigenvalues of DT/^gizo) and are given in (|28ll29p . 

is positive definite, where DPT is computed in the coordinate systems of Ni and 
Ni+i. Notice, to compute DPT we need second order derivatives of the map 
T. We used the C^-Lohner algorithm [WZ2] and the CAPD library [CAPD] in 
order to integrate the second order variational equations for We used the 

Taylor method of the order 20 and a variable time step. No subdivision of the 
sets Ni were necessary, i.e. whole sets Ni were used as an initial condition for 
the routine which computes DPT. The C++ program which verifies the cone 
conditions from this lemma executes within 1 seconds a computer with the Intel 
Xeon 5160, 3GHz processor. I 

7.3 Parameterization of center-unstable and center-stable 
manifolds at cq and Cg, respectively. 

Lemma 22 The following statements hold true. 

• The center-stable manifold of PT at cg can be parameterized as a vertical 
disk in Ng satisfying the cone conditions with respect to the quadratic form 

• The center-unstable manifold of PT at cq can be parameterized as a verti- 
cal disk in Nq satisfying the cone conditions with respect to the quadratic 
form Qq. 

The computer assisted proof of the above lemma is essentially the same as the 
proof of Lemma [17] and Lemma [18] for the Henon map. Therefore we skip 
the details. The C++ program which verifies the assertion executes within 11 
seconds on the Intel Xeon 5160, 3GIIz processor. The most time-consuming 
part is the verification of the existence of covering relations for the projected 
sets in three-dimensional space. 
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We used the Taylor method of the order 20 and a variable time step when 
check the covering relations and when integrate second order variational equa- 
tions. 

Proof of Theorem I19t We conclude the proof as in the Henon map example. 



8 Implementation notes. 

In order to compute bounds for the Henon map and the Poincare map Tp and 
their derivatives we used the interval arithmetic |M], automatic differentiation 
[G] and the C-Lohner algorithm |WZ2] developed at the Jagiellonian University 
by the CAPD group [CAPD] . The C-I--I- source files of the program with an 
instruction how it should be compiled and run are available at [W] . 

The program has been tested under several linux distributions, including 32 
and 64 bits architectures and the gcc compiler versions 4.1.2, 4.2.1 and 4.3.2 on 
the Intel Pentium IV, Intel Core 2 Duo, Intel Xeon and the AMD Quad Core 
processors. 
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